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I review the status of the extraction of NRQCD color-octet J/ j/j production matrix elements. Large theoretical uncertainties in current 
extractions from hadropoduction and photoproduction are emphasized. Leptoproduction of J/ip is calculated within the NRQCD 
factorization formalism. Color-octet contributions dominate the cross section, allowing for a reliable extraction of (Og (*So)) and 

(Og^( 3 Po)). A comparison with preliminary data from the HI collaboration shows that leading order color-octet mechanisms agree 
with the measured total cross section for Q 2 > Am 2 , while the color-singlet model underpredicts the cross section by a factor of 5. 



The NRQCD factorization formalism of Bodwin, 
Braaten, and LepagiJ has emerged as a new paradigm 
for computing the production and decay rates of heavy 
quarkonia. This formalism provides a rigorous theoret- 
ical framework which systematically incorporates rela- 
tivistic corrections and ensures the infrared safety of per- 
turbative calculationsa In the NRQCD factorization for- 
malism, cross sections for the production of a quarkonium 
state H are written as 



a(H) 
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where tuq is the mass of the heavy quark Q. The short- 
distance coefficients, c„, are associated with the produc- 
tion, at distances of order 1/toq or less, of a QQ pair 
with quantum numbers indexed by n (angular momen- 
tum, 2S+1 Lj, and color, 1 or 8). They are computable 
in perturbation theory. In Eq. (0), (0|O^|0) 
matrix elements of NRQCD operators: 



are vacuum 



(0|Of |0> = ]T 5>I4 \H(X) + X){H(X) + X\K n \0), 
x A 

(2) 

where K, n is a bilinear in heavy quark fields which creates 
a QQ pair in an angular-momentum and color configu- 
ration indexed by n. The bilinear combination K) n K, n 
has energy dimension d n . The production matrix ele- 
ments describe the evolution of the QQ pair into a fi- 
nal state containing the quarkonium H plus additional 
hadrons (X) which are soft in the quarkonium rest frame. 
Throughout the remainder of this talk, a shorthand no- 
tation will be used in which the vacuum matrix elements 
are written as (Of l s) { 2S+1 Lj)). 

IThe NRQCD matrix elements obey simple scaling 
laws! with respect to v, the relative velocity of the Q and 
Q in the quarkonium bound state. Therefore, Eq. (Jj]) is a 
double expansion in v and a s . Since the NRQCD matrix 
elements are sensitive only to large distance scales, they 
are independent of the short-distance process in which 
the Q and Q are produced. Thus, the NRQCD matrix 



elements are universal parameters which can be extracted 
from one experiment and used to predict production cross 
sections in other processes. 

Prior to the innovations presented in Ref. 0, most 
J/ip production calculations took into account only the 
hadronization of cc pairs initially produced in a color- 
singlet 3 S\ state, as parameterized by the NRQCD ma- 
trix element (0(^ ( 3 Si)). An important aspect of the 
NRQCD formalism is that, in addition to the color- 
singlet contribution, it allows for the possibility that a 
cc pair produced in a color-octet state can evolve non- 
perturbatively into a. J/ip. The most important color- 
octet matrix elements are (O^ ( 3 S^) , (O^ ( X S )), and 
(Og^( 3 Pj)}, which are suppressed by v 4 relative to the 
leading color-singlet matrix element. They describe the 
non-perturbative evolution of a color-octet cc pair in ei- 
ther a Si, 1 Sq, or 3 Pj angular momentum state into a 
J/^p. Using heavy quark spin symmetry relations, it is 
possible to express all three P-wave matrix elements in 
terms of one: (C 8 ( 3 P 7 )) = (2J + l)(O 8 ( 3 P )} + 0{v 2 ). 
Thus, at this order in the v expansion, there are three 
independent color-octet matrix elements. 

The NRQCD factorization formalism has enjoyed 
considerable phenomenological success. Most notably, 
the large p± production of J/ip at hadron colliders, which 
is underpredicted by over an order of magnitude in the 
color-singlet model (CSM), can be easily accounted for 
by color-octet production mechanismsel The CSM can- 
not account for the observed branching ratios for B —> 
J/ip + X, B — > ip' + X. Inclusiacuof color-octet mech- 
anisms removes this discrepancytffl Color-octet mecha- 
nisms also improve the understanding of J/tp production 
in Z° decaytl 

The NRQCD factorization formalism has yet to be 
conclusively proven as the correct theory of quarkonium 
production. One important test of the formalism is the 
verification of the universality of the NRQCD matrix el- 
ements. Another important prediction, which has yet 
to be verified, is the polarization of J / ip at large p± at 
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hadron colliders. Large p± production of J/ip is domi- 
nated by fragmentation of a (nearly) on-shell gluon into 
a 3 s[ cc pair, which inherits the gluon's transverse po- 
larization. Because of heavy quark spin symmetry, soft 
gluons emitted as the cc pair hadronizes into_the quarko- 
nium state do not dilute this polarizationQcl Therefore, 
at large p±, J/ip are expected to be almost 100% trans- 
versely polarized. This gets significantly Jikited at lower 
due to nonfragmentation productionEjO 
Polarization is a particularly interesting test of 
quarkonium production because it can distinguish be- 
tween NRQCD and the other remaining model of quarka= 
nium production, the color evaporation model (CEM)H3 
In this model, the cross section for producing a J/ip is 
proportional to the total production rate for cc pairs with 
invariant mass less than the DD threshold: 



r2mD j da cE 
vj/-<I> — PJ/ii I dm r .p.- 



drricR ' 
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where pj/^p is a universal factor which represents the frac- 
tion of cc pairs which hadronize into a J/tp. Because 
this model includes color-octet production mechanisms, 
it is consistent_spth hadron collider data as well as data 
from Z° decajej. A prediction of this model is that the 
emission of soft gluons in the hadronization washes out 
any polarization of the cc produced in the short distance 
process. This prediction is at odds with heavy quark 
spin symmetry arguments. The CEM would be a reason- 
able model if the relative velocity, v, in the quarkonium 
state were too large to serve as a useful expansion pa- 
rameter. From the existing estimates of the color-octet 
matrix elements, this does not appear to be the case. 
The leading color-singlet matrix element is well deter- 
mined: (0^( 3 5i)) = 1.1 ± 0.1/, GeV 3 . Extractions 
of (0 8 W ( 3 ^)}, (O 8 j/ *(^o)>, and {O J 8 ^( 3 P )}/m 2 c suffer 
from considerable uncertainty but they are known to be 
of order 10 -2 GeV 3 . These values are consistent with the 
v 4 suppression expected on the basis of NRQCD scaling 
laws. 

To quantitatively predict polarization as a function 
of p±, color-octet matrix elements need to be accu- 
rately determined. At the present time, extractions from 
production at hadron colliders suffer from copsiderable 
theoretical uncertainty. Beneke and Krametll extract 

(Og / ' A ( 3 S'i)) = 1.2jlj$ x 10_2 Gev3 . where statistical er- 
rors, and errors due to variation of the renormalization 
scale and the parton distribution functions, have been 
added in quadrature. The erroLia dominated by scale 
uncertainty. Other extractionljil-a of (0% ( 3 S{)) are 
consistent within these errors. (O^ ^Sq)), (0 8 7/V '( 3 P o )} 
are even more poorly determined since hadroproduction 
is only sensitive to a particular linear combination. The 
extracted value of this linear combination is also much 



more sensitive to the choice of parton distribution £mic=| 
tion and the effects of intital and final state radiatk>nESIi3 
Extractions from lower energy hadroproductiorO suffer 
similar from problems. These calculations are extremely 
sensitive to the charm quark mass because of the small 
energy scale of the process. Corrections due to higher or- 
ders in perturbation theory and higher twist effects are 
expected to be large. 

Recently, progress has been made in analysis of B 
decays an<LZ° decay. A next-to-leading order analysis 
of B decap establishes a bound on the following linear 
combination of color-octet matrix elements: 



(0 8 J/ n i 5o))+3.1 (0 ^ ( 2 3Po)) 



< 2.3 x 10~ 2 GeV 3 .(4) 



The individual matrix elements, (O'^'^Sq)) and 

(Og^( 3 Po)) remain undetermined. Production of J/ip 
in Z° decays is dominated by gluon fragmentation and is 
therefore sensitive to (0%^ ( 3 Si)). A recent analysis!^ 
which sums logarithms of 2Mj i^/Mz and IEji^/Mz to 
all orders, extracts the following effective matrix element: 

(Os^eSi)) = J2(^( 3 S 1 ))Br(H -> J/ip + X) 

H 

— 1.9 ± 0.5 ± 1.0 x 10~ 2 GeV 3 . (5) 

The first error is statistical, the second theoretical. 
The effective (0%^ ( 3 S\)) includes feeddown from higher 
mass charmonium resonances (H) because it is not pos- 
sible to seperate prompt J / ip experimentally. 

Next-to-leading order calculations of photoproduc- 
tion would seem to indicate that NRQCD predicts a large 
excess^ of J/ip at large z which is not observed in the 
dataE3 (Here z — Ej/^/Ey, where energies are measured 
in the proton rest frame.) This led to the conclusion that 
the color-octet matrix elements extracted from the Teva- 
tron were an order of magnitude too large to be consis- 
tent with photoproduction data from the HERA experi- 
ments. However, this enhancement is actually an artifact 
of the next-to-leading order perturbative calculation. It 
is worth examining this process in more detail because 
it is indicative of some of the pitfalls one may encounter 
when trying to compare NRQCD with data. 

The leading order diagrams are shown in Fig. El and 
results in the following perturbative cross section: 



^ = (70*5(1 - z), 
az 

47r 3 a s a 2 
cr = T G{Am z c /s lp ) x 

s fP m c 



(G) 
(7) 



2 



Figure 1: Leading order diagrams for color-octet photoproduction 
of J/rj) 



Figure 2: Next-to-leading order diagrams for color-octet photopro- 
duction of J/ip 



G{x g ) is the gluon distribution function of the proton. 
s-yp is the center of mass energy squared of the photon- 
proton system. The delta function in Eq. (|^) is an artifact 
of the perturbative expansion which will be smeared by 
nonperturbative emission of soft gluons as the cc pair 
hadronizes into the final state hadron. How this non- 
perturbative smearing can be accounted for within the 
NRQCD formalism will be discussed below. For now, we 
concern ourselves with perturbative corrections. 

The next-to-leading order diagrams are shown in Fig. 
||. In the limit in which the final state gluon is soft, the 
cross section is given by: c2l 

da f . ,E|M LO | 2 



dz 



dXgG(Xg)' 



1 QlTX g S'-yp 



xC A g 2 s 



2P-g 



Ami 



P-kg-k {P-kf 



dxG(p/x 



E \Mlo\ 
4m 2 s JP 



■CAOts 



{l~x) 2 {l-z) 



-a - 



Ca& s ln(l - z) 



1-z 



(8) 



Mlo is the leading order matrix element from Fig. 
P, g, and k are the four momentum of the J/i/j, the 
initial state gluon, and the final state gluon, respec- 
tiverly. p — Am^/s^p and x — p/x g . In the last line 
of Eq.(||), we take x = 1 inside the argument of the 
structure function in order to obtain an analytic expres- 
sion. This approximation is valid up to subleading logs 
since the leading logs come from the region x — > z — ► 1. 
In higher orders, we expect to find terms of the form 
a™lri n ~ (1— z) I (1— z) which need to be resummed before 
theory can be compared with experimental data. The di- 
vergence of the cross section as z — > 1 is an artifact of the 
next-to-leading order calculation. Once the resummation 
is performed, the cross section will be well-behaved as 
z -> 1. 

There are also important nonperturbative 
correctionfiJ near 2 = 1. For quarkonium production 
near the boundaries of phase space, it is sometimes the 
case that contributions from NRQCD operators which 



are higher order in v are enhanced by kinematic factors. 
This results in the breakdown of the NRQCD expansion. 
The crux of the problem is that in the perturbative QCD 
part of the matching calculation one uses twice the heavy 
quark mass instead of the quarkonium mass to compute 
the phase space for the production of the quarkonium 
meson. The difference between 2mg and Mh is a v 2 cor- 
rection, which is ignored in leading-order calculations. 
However, at the boundaries of phase space this differ- 
ence becomes important, and it is necessary to sum an 
infinite number of NRQCD matrix elements. This resum- 
mation leads to a universal distribution function called a 
shape function which replaces the delta function in Eq. 
(||). Because of the universality of the shape functions, it 
may be possible in the future to test NRQCD by compar- 
ing shape functions extracted from different quarkonium 
production processes. This would require more precise 
data than is currently available. 

As z — > 1, the photoproduction cross section becomes 
sensitive to the effects of soft gluon radiation. Within 
pertubation theory this sensitivity is signalled by the 
presence of large logarithms which must be resummed. 
The v expansion in Eq. (|l|) also breaks down and an in- 
finite set of operators become relevant. For this reason, 
it is impossible to extract NRQCD matrix elements from 
this distribution. When leptoproduction is discussed be- 
low, the analysis will focus on observables for which shape 
function corrections are not kinematically enhanced, be- 
cause it is only from these distributions where one can 
hope to obtain a reliable extraction of NRQCD matrix 
elements. 

In RefEi, leptoproduction of J/tp is examined and 
found to be a useful process from which to extract 
NRQCD matrix elements. This process is similar to pho- 
toproduction, however, now the photon is off-shell. At 
leading order, the cross section is again given by the di- 
agrams in Fig. |l|: 



da 
dQ 2 



= dy dx g G(x g ) 5(x g ys - (2m c ) 2 - Q 2 ) 
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2a s {y)c 



1 + (1 - yf 



Q 2 (Q 2 + (2m c 



3Q 2 + 7(2m c ) 2 (OtfPo)} 



< 2 + (2m c ) 2 
8(2m c ) 2 Q 2 



(^( 3 ^o)) 



x g s(Q 2 + (2m c ) 2 ) m 2 



(9) 



where s is the electron-proton center-of-mass energy 
squared. The momentum fraction of the virtual pho- 
ton relative to the incoming lepton is y = P p ■ q/P p ■ k, 
where P p is the proton four-momentum, q is the pho- 
ton four-momentum, and k is the incoming lepton four- 
momentum, and Q 2 = —q 2 . 

There are several advantages to extracting NRQCD 
matrix elements from this process. First of all, note 
that the relative importance of (O^( 1 5 )> and (Of( 3 P )) 
changes as a function of Q 2 . Thus, it is possible to fit 
the differential cross section as a function of Q 2 and ex- 
tract both of these matrix elements seperately. Second, 
the scale for the coupling is set by /i 2 w Q 2 + (2m c ) 2 , 
so as Q 2 increases, perturbative corrections should be 
increasingly suppressed relative to photoproduction. In 
RcfE3, next-to-leading order graphs with real gluon emis- 
sion were computed and found to give a small contribu- 
tion to the total cross section for Q 2 > 4 GeV 2 . The 
uncertainty due to varying the renormalization scale is 
less than 10% for Q 2 > 4 GeV 2 . This is a great deal 
smaller than the corresponding uncertainty in hadropro- 
duction. This indicates that higher order perturbative 
corrections to the leading order color-octet production 
mechanism will be small. Higher twist corrections to the 
parton model are also expected to be small as they are 
suppressed by powers of Q 2 . 

The leading color-singlet mechanism requires the 
emission of an additional hard gluon. (See the first di- 
agram in Fig. 0.) This a s (Q) suppression compensates 
for the u 4 suppression of the color-octet matrix elements, 
and it turns out that for Q 2 > 4 GeV 2 , the color-octet 
contribution is roughly 4 times bigger than the leading 
color-singlet process. However, color-singlet mechanisms 
will dominate in the inelastic (z < 1) region, where there 
is no 0(a s ) color-octet contribution. In the elastic re- 
gion, one also expects production of J/ip via diffractive 
processes. At large Q 2 , diffractive leptoproduction can 
be studied using perturbative QCD. In Ref. E3, a pertur- 
bative analysis of diffractive leptoproduction predicts the 
cross section to fall as 1/(Q 2 + (2m c ) 2 ) 3 , as compared to 
1/(Q 2 + (2m c ) 2 ) 2 for the color-octet mechanism. There- 
fore, at sufficiently large Q 2 , the diffractive contribution 
should be negligible correction to our calculation. 

One must also consider the possibility of kinematic 
enhancement of higher order w 2 _corrections. An analysis 
of the Q 2 distribution in Ref. E3 shows that the higher 



order in v 2 corrections associated with the shape func- 
tion are suppressed in the large Q 2 limit. Therefore, this 
distribution is calculable in NRQCD. However, the z dis- 
tribution still suffers from large endpoint corrections, and 
a perturbative calculation of this distribution cannot be 
compared with experiment near z = 1, even at large Q 2 . 

The largest errors in the calculation of the Q 2 dis- 
tribution are associated with uncertainty in the charm 
quark mass. Varying m c between 1.3 GeV and 1.7 GeV, 
results in an error of ^25% a ^ Q 2 = ^ GeV 2 . This error 
decreases slightly as Q 2 is increased. 

Once color-octet matrix elements are extracted, the 
polarization of the J ftp can be predicted without intro- 
ducing any new parameters. Calculations of the polar- 
ization of J I ip produced via the leading order color.=qctet 
mechanism as a function of Q 2 are given in Ref. E3. A 
precise measurement of the polarization of leptoproduced 
J ftp could provide an excellent test of the NRQCD fac- 
torization formalism. 

Preliminary results from HERApJ indicate that 
da/dQ 2 is well described by NRQCD for Q 2 > (2m c ) 2 . 
This is in contrast with CSM calculations which undcr- 
predict the cross section by a factor of ~ 5. So this ap- 
pears to be strong evidence in favor of color-octet mech- 
anisms as understood within the NRQCD factorization 
formalism. However, the shape of the rapidity distribu- 
tion at large rapidity is not well reproduced by NRQCD. 
This could be a consequence of the effects of soft gluon 
emission which are not accounted for in the leading order 
calculations of Ref. El. Analysis of soft gluon effects and 
extraction of NRQCD matrix elements will be performed 
in a future publications 
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